Measurement of the flux of atmospheric muons with the CAPRICE94 apparatus 
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A new measurement of the momentum spectra of both pos- 
itive and negative muons as function of atmospheric depth was 
made by the balloon-borne experiment CAPRICE94. The 
data were collected during ground runs in Lynn Lake on the 
19-20th of July 1994 and during the balloon flight on the 
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8-9th of August 1994. We present results that cover the mo- 
mentum intervals 0.3-40 GeV/c for /i~ and 0.3-2 GeV/c for 
/i+ , for atmospheric depths from 3.3 to 1000 g/cm^, respec- 
tively. Good agreement is found with previous measurements 
for high momenta, while at momenta below 1 GeV/c we find 
latitude dependent geomagnetic effects. 

These measurements are important cross-checks for the 
simulations carried out to calculate the atmospheric neutrino 
fluxes and to understand the observed atmospheric neutrino 
anomaly. 
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I. INTRODUCTION 

Recent atmospheric neutrino observations by Super- 
Kamiokandc and, with lower statistics, Soudan-2 and 
MACRO cohaborations have been used as evidence of 
neutrino oscillation. The observed rate of neutrino inter- 
actions were compared to the rates calculated using the 
neutrino fluxes derived from atmospheric cascade simu- 
lations It was found that the observed number of 
events induced by muon neutrinos is too few compared 
to that from simulation. While several aspects of the 
measurements, such as the zenith arigular dependence of 
the observed atmospheric neutrinos ^ , strongly point to- 
ward the hypothesis of neutrino oscillations, the precise 
determination of the allowed and excluded regions in the 
oscillation parameter space rely heavily on the compari- 
son of measurements with calculations . 

Recently Gaisser et al. |^ compared different calcula- 
tions of the atmospheric neutrino flux and concluded that 
the differences can be attributed to three main effects. 
First, the different parametrizations used of the particle 
production and pion yield in hadronic interactions of the 
primary cosmic rays with air nuclei. Second, the abso- 
lute energy spectra of the primary cosmic rays (protons 
and helium nuclei). Third, the solar modulation and ge- 
omagnetic effects. These differences reflect the existing 
experimental uncertainties on these topics. The avail- 
able accelerator data on the pion yield in interactions of 
hadrons with air nuclei are limited and the results used in 
the simulations do not cover the whole phase space par- 
ticularly at low values of the Feynman x. Older measure- 
ments of the primary cosmic ray spectra differ by as 
much as 40% from more recent observations |^-[l0[] which 
are in agreement at the level of 10%-20%, compatible 
with the statistical and systematic uncertainties of the 
measurements. Besides the uncertainties on the absolute 
values, the primary cosmic ray fluxes, particularly be- 
low 10 GeV/n, are affected by the periodic solar activity 
and the geomagnetic field. Furthermore, the geomagnetic 
field affects the distribution of the cascade in the atmo- 
sphere and these effects are not accounted for in most 
of the simulations which use a 1-dimensional approach, 
that is all secondaries are assumed to be collinear with 
the direction of the primary particles. Gaisser et al. |^] 
concluded that the calculated neutrino interaction rates 
are uncertain at the level of ±30%. They also concluded 
that because of the cancellation of errors, the ratio of i'^ 
to Ve has a smaller uncertainty, on the order of ±5%. 

The majority of the Super-Kamiokande events are in 
the Sub-GeV region [Q where the above uncertainties 
are the most critical. Measurements of the flux of atmo- 
spheric muons provide a powerful tool for cross-checking 
the cascade simulations. Nearly all the sub-GeV neutri- 
nos events originate from the tt — > — > e decay chain. 
Neutrinos with an energy of 1 GeV are generated, on av- 
erage, by muons with an energy of about 3 GeV. This 
interaction typically takes place at altitudes between 12 



and 26 km (200 to 20 g/cm^). Hence, muon measure- 
ments must be performed over a broad energy range, 
from a few hundred MeV to tens of GcV, and over an 
extended range of atmospheric depths in order to be suit- 
able for use in cross-checking. 

We present in this paper results on the muon spec- 
tra in the atmosphere obtained with the CAPRICE94 
instrument from ground (360 m) to float (36-38 km) al- 
titude. First results on this analysis were reported ear- 
lier jl|jl^. In ||ll| we compared our measured fluxes 
to results of simulations and concluded that the calcula- 
tions overestimated our measured muon fluxes; the differ- 
ences depending on momentum and atmospheric depth. 
Theoreticians are introducing changes in the simulations 
procedures, e.g. changing from one-dimensional to three 
dimensional interaction models, to account for these dis- 
crepancies. Therefore, in this paper we do not make any 
comparisons with simulation results. In this paper we 
provide details of the data analysis and present our final 
results. We describe the detector system in this experi- 
ment in section 2, the data analysis in section 3 and the 
results and discussion comprise section 4. 

It is worth mentioning that the primary cosmic ray hy- 
drogen and helium spectra also were measured with 
the CAPRICE94 instrument. These can be used as the 
input spectra for the cascade simulations in order to re- 
duce the overall systematic uncertainties associated with 
the comparison of observed and calculated muon fluxes. 



II. DETECTOR SYSTEM 

Figure I shows the NMSU-WiZard/CAPRICE spec- 
trometer that was operated on ground at Lynn Lake, 
Manitoba, Canada (56.5° North Latitude, 101.0° West 
Longitude, 360 m altitude), in summer 1994. The pay- 
load was flown by balloon from Lynn Lake to Peace River, 
Alberta, Canada (56.15° North Latitude, 117.2° West 
Longitude), on August 8-9, 1994. The balloon reached 
a altitude of 36.0 km, corresponding to an atmospheric 
pressure of 4.5 mbar, in about three hours of ascent and it 
floated at altitudes ranging from 38.1 to 36 km, i.e. resid- 
ual atmosphere of 3.3-4.6 g/cm^, for about 23 hours. The 
apparatus included from top to bottom: a Ring Imaging 
Cherenkov (RICH) detector, a time-of-flight (ToF) sys- 
tem, a superconducting magnet spectrometer equipped 
with multiwire proportional chambers (MWPC) and drift 
chambers (DC) and a silicon-tungsten imaging calorime- 
ter. 

The 51.2x51.2 cm2 RICH detector |l|l, used a 1 cm 
thick solid sodium flouride (NaF) radiator with a thresh- 
old Lorentz factor of 1.5, and a photosensitive MWPC 
with pad readout to detect the Cherenkov light image 
and hence measure the velocity of the particles. 

The time-of-flight system had two layers one above and 
one below the tracking stack, each layer made of two 1 cm 
thick 25 X 50 cm^ paddles of plastic scintillator. Each pad- 



2 



die had two 5 cm diameter photomultiplier tubes, one at 
each end. The distance between the two scintillator lay- 
ers was 1.1m. The time-of- flight system was used to give 
a trigger, to measure the time-of-flight and the ionization 
(dE/dX) losses of the particles. The trigger was a four- 
fold coincidence between two photomultiplier tubes with 
signals in the top paddle and two in the bottom scintil- 
lator paddle. The threshold of each photomultiplier tube 
was set high enough to eliminate noise, but low enough to 
provide an efficiency of nearly 100% to trigger minimum 
ionizing particles. 

The spectrometer consisted of a superconducting mag- 
net and a tracking device equipped with multiwire pro- 
portional chambers and drift chambers . The magnet 
consisted of a single coil of 11 161 turns of copper-clad 
NbTi wire. The outer diameter of the coil was 61 cm 
and the operating current was 120 A, producing an in- 
homogeneous field of about 4 T at the center of the coil. 
The spectrometer provided 19 position measurements (12 
DC and 7 MWPC) in the direction of maximum bending 
(x) and 12 measurements (8 DC and 4 MWPC) along 
the perpendicular direction (y). Using the position infor- 
mation together with the map of the magnetic field, the 
rigidity of the particle was determined. 

Finally, the electromagnetic calorimeter consisted 
of eight 48x48 cm^ planes of silicon strip (3.6 mm wide) 
detectors with both x and y readouts. These silicon 
planes were interleaved with layers of tungsten convert- 
ers, each one radiation length thick. Taking into ac- 
count all the material, the calorimeter had a total thick- 
ness of 7.2 radiation lengths and 0.33 nuclear interaction 
lengths. The calorimeter provided topological informa- 
tion on both the longitudinal and lateral profiles of the 
particle's interaction as well as a measure of the total 
energy deposited in the calorimeter. 



III. DATA ANALYSIS 

The analysis was based on 85800 seconds of data taken 
at ground, 10000 seconds during the ascent of the payload 
and 60520 seconds of data taken at float altitude under 
an average residual atmosphere of 3.9 g/cm^. 



of the RICH to reject proton events. At increasing atmo- 
spheric depths, the abundance of the proton component 
decreased to a few percent of the positive muon compo- 
nent at ground level. For negatively charged particles 
at float, the electron component was the dominant back- 
ground. The electron background rapidly decreased with 
increasing atmospheric depth becoming smaller than the 
muon component by 200 g/cm^ and a small fraction of 
the muon component at ground level |l^,|l^ • Because of 
this varying background different selection criteria were 
used for /i+ and /i^ and for ground, ascent and float 
data to maximize the efficiency while keeping the reject- 
ing power for background events at an appropriate level. 

The acceptance of the instrument allowed for muons 
with a range of zenith angles to be measured. The maxi- 
mum angle was 20 degrees and the mean of the distribu- 
tion was at 9 degrees. Figure |2| shows the cosine zenith 
angle distribution of for muons of both signs selected be- 
tween 0.15 and 2 GV/c at ground level, during the ascent 
and at float. The distributions have been normalized for 
the total number of events. No significant change is found 
in the zenithal angle distribution. It is worth pointing out 
that the distribution narrows as the rigidity increases. 

The selections used for identifying muons in the 
ground, ascent and flight data are summarized in Ta- 
ble I and described below. The figures describing several 
of these selections show float data since at float the back- 
ground of other particles was the largest. 



1. Tracking 

The tracking information was used to determine the 
rigidity of the particles. In this work the trajectory was 
determined by fitting only the information from the drift 
chambers. This made it possible to use the MWPC sys- 
tem for the efficiency estimation of the drift chambers 
To achieve a reliable estimation of 



(see section III C 1 



the rigidity, a set of conditions were imposed on the fitted 
tracks: 

1. At least 9 (out of 12) position measurements in the 
X direction and 5 (out of 8) in the y direction were 
used in the fit. 



A. Particle selection 



2. There should be an acceptable chi-square for the 
fitted track. 



The CAPRICE94 instrument was well suited to mea- 
sure the muon spectra and charge ratio in the atmo- 
sphere against a background of electrons, protons and 
heavier particles. The background in the muon sample 
depended strongly on the atmospheric depth. At float 
the dominant background for positive particles was pro- 
tons, which outnumbered the positive muons by about 
a factor of 1000. The upper end of the energy interval 
for the fi^ measurement was determined by the ability 



3. The estimated error on the deflection should be < 
0.04 (GV/c)-i. 

These conditions also eliminated events with more than 
one track in the spectrometer. 



2. Scintillators and time-of-flight 

The ionization (dE/dX) loss in the top ToF scintilla- 
tor was used to select minimum ionizing singly charged 
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particles by requiring a measured dE/dX of less than 1.8 
times the most probable energy loss by a minimum ion- 
izing particle (mip). 

Downward moving particles were selected using the 
time-of-flight information. The time-of-flight resolution 
of 280 ps, which was small compared to the flight-time 
of more than 4 ns, assured that no contamination from 
albedo particles remained in the selected sample. The 
particle's velocity {(3) reconstructed from the time-of- 
flight information was used to select muons against a 
background of pions, protons and heavier particles. Fig- 
ure shows P obtained from the time-of-flight informa- 
tion for positive particles collected at float as a function of 
rigidity. Muons were selected as particles above the solid 
line. The ToF (3 selection was used for identification 
for the ascent and float portion of the data and for the 
/i~ selection during float. 



3. The calorimeter. 

With a depth of 7 radiation lengths the calorimeter 
could identify non-interacting particles in a background 
of electrons above a momentum of a few hundred MeV/c. 
The selection was performed by requiring that the topo- 
logical information of the signal was consistent with that 
of a single track. This was accomplished by imposing 
an upper limit to the number of hits along the track in 
the calorimeter. This selection reached its highest ef- 
ficiency above 1 GeV/c where electromagnetic showers 
were well defined in the calorimeter. Below 1 GeV/c a 
non-negligible electron contamination was present. To 
further reduce this background, another selection crite- 
rion, based on the total detected energy in the calorime- 
ter divided by the momentum, was used. An upper limit 
for this quantity equal to 60 mip/(GeV/c) was applied. 
Figure ^ shows this quantity as a function of rigidity 
for the float data. The two dense band are due to non- 
interacting particles. Recall that muons above a few hun- 
dred MeV/c are minimum ionizing particles that release 
a nearly constant amount of energy in the calorimeter. 
Below 300 MeV/c the calorimeter selection was not used 
because of the low efficiency. At ground level e^ amount 
to less than 0.1% of the muon component above 3 GV/c 
p6| , p7t and, consequently, the calorimeter muon selection 
criteria were not used above this rigidity. 

The calorimeter selection criteria also rejected particles 
that interacted in the calorimeter, namely pions, protons 
and heavier nuclei, hence contributing to reduce their 
contamination in the muon sample. However, the re- 
jection factor for these particles was small because the 
calorimeter depth was only one third of a nuclear inter- 
action length. Therefore, hadrons were removed using 
selection criteria based on time-of-flight, scintillator and 
RICH information. 



4. The RICH 

The RICH was used to measure the Cherenkov angle 
of the particle and thereby its velocity. The Cherenkov 
angle was reconstructed from the geometrical distribu- 
tion of the signals in the pad plane (for a description of 
the reconstruction method see ||l^). To correctly use the 
RICH information, a set of conditions was applied on the 
RICH data They were: 

1. Ionization from charged particles produced signif- 
icantly higher signals than converting Cherenkov 
photons. To reject events with multiple charged 
particles traversing the RICH, we required that an 
event contained only one cluster of pads with high 
signals. 

2. A good agreement between the particle impact po- 
sition as determined by the RICH and the track- 
ing system was required. The difference in x and 
y should be less than three standard deviations, 
which was rigidity dependent but typically less 
than 5 mm. 

3. More than 3.5 pads with signals from converted 
Cherenkov photons were required in the fit for the 
ground data and more that 7.5 for the ascent and 
float data. 

4. The reconstructed Cherenkov angle should not de- 
viate by more than three standard deviations from 
the expected Cherenkov angle for muons. 

Criterion 1 was introduced to eliminate events with more 
than one charged particle crossing the RICH and this 
condition was applied over the entire data sets both for 
/i+ and fi^ . Criterion 2 eliminated events that scattered 
in the RICH electronics. Criteria 3 and 4 were used to 
separate muons from the other particles. Figure ^ shows 
the measured Cherenkov angle for flight data events se- 
lected with RICH criteria 1, 2 and 3. The bands cor- 
responding to the different particles are clearly visible. 
The solid lines indicate the muon selection based on the 
Cherenkov angle. However, it is important to point out 
that the solid lines in the figure are only indicative of 
the selection, since the RICH selection was done on an 
event by event basis. For each event the Cherenkov angle 
was obtained and the resolution of this depended on the 
incident angle of the particle . 

These four criteria were used for selecting muons of 
both signs at float and for selecting /i+ at ground and 
during the ascent. Ground and ascent /i~ were selected 
differently: the selection was based upon the dE/dX and 
calorimeter criteria except below 70 g/cm^ of residual 
atmosphere where contamination of int eraction products 
from the payload structure (see section HI B 4 ) was non- 
negligible. For these small depths the RICH criterion 1 
also was used. In addition the full RICH selection with 
criteria 1 to 4 was used below 500 MeV/c since, at low 
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velocities, the RICH was able to separate muons from 
pions and electrons. 

5. The bar 

A 17 kg, 1.2 m long aluminum bar with a 7 kg steel 
hook in the centre was used to connect the payload to 
the balloon during the flight. This bar was situated 2.3 
m above the RICH. Hadrons had a non-negligible proba- 
bility to interact with the material of the bar and produce 
secondary particles that could be detected as muons in 
the apparatus. Hence we chose to reject all particles with 
extrapolated trajectories that crossed the bar. This pro- 
cedure resulted in a reduction of the geometrical factor 
by about 10% as can be seen in Figure |[ 

B. Background rejection 

The various sources of background in the muon analy- 
sis are described below along with the rejections criteria 
and surviving contamination. This information is sum- 
marized in Table O. 



1. Electron background 

The calorimeter muon selection gave an electron rejec- 
tion factor that increased from 30 at 0.3 GeV/c to more 
than 1000 above 1 GeV/c. The RICH separated muons 
from electrons with an electron rejection factor of more 
than 100 at 0.1 GeV/c decreasing to 10 at 0.3 GeV/c 
and to 1 at 0.7 GeV/c. Since the electron flux is higher 
than the muon one only at high altitudes where the elec- 
tron primary component dominates and is about a fac- 
tor four larger than the muon flux |l7|j2C|| , the electron 
and positron contaminations surviving the muon selec- 
tion were assumed negligible. 

2. Proton background. 

The ToF (3 rejection factor for protons below 1 GeV/c 
was greater than 4000 decreasing to about 300 at 
1.2 GeV/c. The RICH proton rejection factor was greater 
than 2000 below 1.5 GeV/c, about 1000 at 2 GeV/c and 
decreased to about 1 at 5 GeV/c. Below 1.5 GeV/c, 
the combined effect of the time-of-flight and RICH selec- 
tion criteria reduced the proton contamination to a neg- 
ligible fraction of the selected muon sample. At higher 
momenta, because of the strong variation of the pro- 
ton flux, the contamination was dependent on the atmo- 
spheric depth. At float altitude the proton contamination 
became increasingly important above 1.5 GeV/c and it 
dominated the muon sample above 2.5 GeV/c. Hence, 



we conservatively limited the positive muon measure- 
ments to a momentum range between 0.15 and 2 GeV/c. 
In this range the proton contamination was negligible 
at all atmospheric depths except at float. In the float 
data we assumed, as a worst case, that all the singly 
charged particles were protons. Applying the rejecting 
power of the calorimeter, dE/dX, ToF beta and RICH 
to this sample we found that a small (less than 20%) 
proton contamination survived the muon selection in the 
rigidity range 1.5 to 2 GeV/c. This contamination was 
subtracted from the positive muon sample at float. For 
ground data the range for /i+ was extended to higher mo- 
menta as presented in p^ , p^ . The proton contamination 
was negligible at ground level below 3 GeV/c, while above 
3 GeV/c it was calculated and subtracted from the posi- 
tive muon sample. This calculation was made by rescal- 
ing the number of interacting particles in the calorimeter 
with factors obtained from data at float. The proton can- 
didates were selected if they had a hadronic interaction in 
the calorimeter. The contamination from muons in the 
interacting proton sample was studied using negatively 
charged ground data, while the contamination from elec- 
tromagnetic showers was negligible at momenta greater 
than 3 GeV/c p6| , p7t . The efficiency of this proton se- 
lection was estimated using a sample of singly charged 
particles at float that was assumed to be composed of 
protons. 

In summary, during the ascent positive muons were se- 
lected free of proton contamination up to 2 GeV/c. At 
float a small proton contamination was subtracted from 
the positive muon sample between 1.5 and 2 GeV/c while 
at ground the range was extended to much higher mo- 
menta subtracting the small estimated proton contami- 
nation 

3. Heavier nuclei background. 

The case of deuteron background was very similar to 
the proton one. Helium and heavier nuclei were mainly 
rejected by the dE/dX selection. The remaining fraction 
was eliminated using the other selection methods and 
their contamination was determined to be negligible. 

4. Meson background. 

Because mesons (pions and kaons) produced in the at- 
mosphere above the payload decay rapidly, they repre- 
sented a small (< 2%) fraction of events compared to the 
muon flux (e.g. see [^). However, pions produced in the 
payload could still present a non-trivial contamination. 
We have undertaken a careful analysis of the local pion 
background in order to quantify their abundance. 

Evidence for a pion contamination is visible at very 
low rigidities in Figure |. These pions, because of their 
low energies, were presumably produced locally in the 
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RICH or in the dome (part of the gondola above the 
RICH). However, it is important to stress that no RICH 
conditions were apphed to select the events in Figure ||: 
the RICH selection rejected multiply charged tracks, par- 
ticles interacting in the RICH as well as pions below 
500 MV/c. Also the dE/dX selection rejected multiply 
charged tracks and this leaves the dome as the source of 
pions produced in the interaction with high energy cos- 
mic rays. The interactions occurred at such an angle that 
the high energy secondaries missed the active volume of 
the instrument, but a low energy pion was produced at 
a large angle and passed through the instrument appear- 
ing like a singly charged particle rather than part of a 
shower. In studying the pions below < 0.2 GV/c, where 
they can be identified using the time-of- flight, we found 
no evidence of a preferred incoming direction. This is 
consistent with our interpretation of these events since 
the distribution of material in the gondola was symmet- 
ric with respect to the azimuthal angle. The locally pro- 
duced pion flux entering the apparatus should decrease 
quickly with energy because of the emission towards the 
forward direction and the lower probability of having only 
one charged particle with high momentum traversing the 
detectors. 

To test the correctness of this conclusion the following 
approach was adopted. Events were selected from float 
data with: 

1. multiple pad ionization clusters in the RICH; 

2. multiply charged signal in the top scintillator; 

3. rigidity (R) interval: 0.5 < i? < 2 GV/c; 

4. Cherenkov angle of a /? « 1 particle; 

5. time-of-flight of a /3 w 1 particle; 

6. dE/dX signal of a minimum ionizing particle in the 
bottom scintillator. 

With these criteria /? ~ 1 particles belonging to a shower 
initiated closely to the top of the apparatus were selected, 
hence with high probability secondary pions, which en- 
tered the calorimeter as singly charged. The criteria 3 to 
6 ensured that the selected events were indeed /3 « 1 par- 
ticles. This sample of locally produced pions was used to 
estimate their interaction probability in the calorimeter 
(or, more precisely, the probability of selecting interact- 
ing pions in the calorimeter). Then, the same hadronic 
interaction selection for the calorimeter was applied to 
the muon sample selected without using the calorimeter 
muon criteria. Two samples were obtained, one for posi- 
tive and one for negative pions. These events (about 100 
in total) were visually scanned with a graphic program. 
In this way misidentified muons and, especially, electrons 
were rejected from the samples. Then, the remaining 
pion numbers were rescaled by the calorimeter selection 
efficiency thus obtaining the number of pions in the muon 
samples. The result was that at float altitude pions could 



account for a maximum of 20% between 0.5 and 1 GeV/c 
and for less than 10% above 1 GeV/c of the muon flux, 
irrespective of the sign of the charge. We give this re- 
sult as an upper limit because the procedure is likely to 
overestimate the number of pions (some of the selected 
pions could be muons, etc.). For this reason, it was not 
subtracted from the muon flux and it should be consid- 
ered as a systematic uncertainty of the measurement. It 
is important to stress that this uncertainty, due to the 
similar pion contamination in the /x"*" and ii~ samples, 
affects the /i+ to /i~ ratio less than the corresponding 
fluxes. 



5. Conclusion on background 

Clean /i+ and fi^ samples were selected from 0.15 to 
0.4 GeV/c. Above this momentum a non-negligible con- 
tamination of locally produced pions could be present. 
For the float data, this contamination was less than 20% 
of the muon flux between 0.5 and 1 GeV/c decreasing to 
less than 10% above 1 GeV/c. For larger atmospheric 
depths the locally produced pion flux decreased quickly 
due to the decrease of the interacting proton and helium 
nuclei, specially at large zenith angle, while the muon flux 
increased with increasing depth at least up to 100 g/cm^. 
Hence the locally produced pion contamination was as- 
sumed negligible at all depths except at float. 

At float a small proton contamination was subtracted 
from the positive muon sample between 1.5 and 2 GeV/c. 
At ground the momentum range was extended to much 
higher value with a subtraction of a small proton con- 
tamination. 



C. Efficiency determination 

In order to accurately determine the fluxes of the vari- 
ous types of particles, the efficiency of each detector was 
carefully studied using both ground and flight data. To 
determine the efficiency of a given detector, a data set 
of muons was selected by the remaining detectors. The 
number of muons correctly identified by the detector un- 
der test divided by the number of events in the data set 
provided a measure of the efficiency. This procedure was 
repeated for each detector. The efficiency of each detec- 
tor was determined as a function of rigidity in a number 
of discrete bins and, then, parameterized to allow an in- 
terpolation between bins. This parameterization intro- 
duced a systematic error on the efficiency of each detec- 
tor. Since the parameters were correlated, the error on 
the efficiency was obtained using the error matrix of the 
fit for each detector when correcting the measured ffux 
for the detector efficiencies. 
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1. Tracking efficiency 

The drift chamber tracking efRciency was obtained us- 
ing negative singly charged particles selected by the other 
detectors similarly as done in |Q. A sample of singly 
charged particles was selected by requiring a single ion- 
ization cluster in the RICH and a dE/dX signal in the top 
scintillator typical of a minimum ionizing particle. From 
this sample, negatively charged events were selected by 
requiring a negative deflection from the fit to the MWPC 
trajectory measurements. The contamination of spillover 
protons was eliminated by requiring that the measured 
impact positions in the RICH and calorimeter agreed 
with the positions as obtained by extrapolating the par- 
ticle trajectory derived from the MWPC fit. The result- 
ing sample of negative singly charged particles was used 
to determine the efhciency of fitting tracks in the drift 
chamber system. The solid line in Figure |^ shows the 
tracking efficiency at float altitude. The efficiency varied 
from ground to float altitude. At ground it was ~ 95% 
above 1 GeV/c, then just after the launch it was ~ 87% 
increasing to ~ 93% at float altitude. 

Biases in the efficiency sample were studied using pro- 
tons (see [Q). It was found that the criteria used for 
fitting the tracks using the MWPC slightly reduced the 
number of scattered tracks in the sample. In order to ac- 
count for this reduction a systematic uncertainty of 2% 
was introduced JTzt . 

Possible charge sign dependence of the efficiency was 
studied using both the flight data and the data taken on 
the ground before the flight. No significant dependence 
was found above 0.3 GV/c 0. 

2. Scintillator ejficiency 

The dE/dX and (3 selections were studied using neg- 
ative events with a minimum ionizing pattern in the 
calorimeter. A clean RICH signal also was required to 
reject interaction products from the sample. The dot- 
ted lines in Figure ^ show the scintillator efficiency in its 
two selections. The efficiency was studied using ground, 
ascent and fioat data and no variation was found. 



3. Calorimeter ejficiency 



since the calorimeter performances were stable over a pe- 
riod of months. 



4- RICH efficiency 

On the ground and in the first part of the ascent the 
RICH efficiency was obtained by selecting negative singly 
charged particles, which did not interact in the calorime- 
ter. At float altitude a large background of interaction 
products did not permit us to select an unbiased clean 
sample of muons, hence the efficiency obtained from the 
ground data was used. This procedure was validated by 
comparing the RICH efficiency for selecting electrons at 
ground and at float since an unbiased clean sample of 
electrons could be selected using the calorimeter. The 
RICH electron selection criteria is similar to that for 
muons, namely, same requirements on the impact po- 
sition, number of effective pads and Cherenkov angle but 
using the theoretical electron Cherenkov angle. It was 
found that the RICH electron efficiency for float data 
reproduced the electron efficiency of ground electrons in- 
side an uncertainty of about 5%. Therefore, it was rea- 
sonable to make use of the muon RICH efficiency as ob- 
tained from the ground data (dashed-dotted line in Fig- 
ure 0) for the flight data. 

D. Geometrical factor 

The geometrical factor, determined with simulation 
techniques pl| , is shown in Figure ^ for ground and flight 
data. The difference at low deflections (high rigidities) 
for the two sets of data is due to the additional geomet- 
rical constrain imposed due to the bar. 

The geometrical factor was cross checked with two 
other methods. One adopted the same approach as pre- 
sented in pl| ] using, however, a different method to trace 
the particles: the track fltting algorithm used in the anal- 
ysis also was used to trace the particle through the spec- 
trometer. This method gave the same results within 1%, 
at all rigidities. The second used a numerical integration 
calculation of the geometrical factor that agreed with the 
previous results within 2% above 0.5 GV/c and within 5% 
below 0.5 GV/c. 



The calorimeter selection efficiency, shown as a dashed 
line in Figure ^, was obtained using ground data. The re- 
sult was cross checked with a simulation of the calorime- 
ter and an excellent agreement was found. The 
calorimeter efficiency also was studied with flight data. 
The calorimeter only is able to separate muons from elec- 
trons above 0.5 GeV/c, hence the presence of a larger 
electron contamination had to be taken into account. In- 
side the errors a good agreement was found as expected 



E. Systematic uncertainties 

Systematic uncertainties originating from the determi- 
nation of the detector efficiencies were included in the 



tables and data points as discussed in session HI C . An- 
other possible systematic error was related to the effi- 
ciency of the trigger system. The fraction of each trigger 
combination was compared with the simulated fraction 
taking into account the position of each paddle and the 
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magnetic field. The excellent agreement between the sim- 
ulated and experimental fractions permitted us to con- 
clude that a possible systematic error due to a geometri- 
cal inefficiency of the trigger was less than 1%. 

The residual atmosphere above the gondola was mea- 
sured by two pressure sensors owned and calibrated by 
the National Scientific Balloon Facility. The two mea- 
surements did not coincide: their difference increased 
with altitude from less than 1% to about 10% at float. 
We interpreted this difference as the systematic uncer- 
tainty on the atmospheric depth. This uncertainty does 
not affect the measurement but has to be taken into ac- 
count when comparing the measured spectra with the 
simulated ones. 



From the discussion in section HID we conclude that 



the systematic error due to the geometrical factor calcu- 
lation was less than 5% between 0.3 and 0.5 GV/c and 
less than 2% for rigidities higher than 0.5 GV/c. For the 
geometrical factor calculations it was assumed that there 
was no variation of the muon intensity over the accep- 
tance angle. The effect on the geometrical factor due to 
the intensity variation was examined using the measured 
muon spectra and the observed zenithal distribution 
in our apparatus in the rigidity range 0.2 to 1.5 GV/c at 
ground. We found that the calculated geometrical factor 
would be reduced by about 3%. 

We decided to assign a systematic uncertainty of 5% to 
the RICH muon efficiency at fioat to account for possible 
variations in the RICH performance between ground and 
fioat. 

The tracking muon selection efficiency varied with time 
during the ascent. We determined the efficiency for seven 
time bins and we found that they could be grouped in 
three in which the efficiency could be assumed constant. 
Since the efficiency above 1 GV/c varied from about 87% 
at launch to 93% at float we believe that the systematic 
uncertainty of this procedure is less than 6%. 

Assuming that the systematic uncertainties discussed 
above are uncorrelated, we estimated an overall sys- 
tematic uncertainty, which is momentum dependent, for 
ground muons decreasing from ~ 6% at 0.3 GeV/c to 
about 2% above 1 GeV/c. It is worth pointing out that 
the ground muon fluxes measured by the CAPRICE94 
apparatus agree at the level of 3% with the measure- 
ments from the C APRICE97 experiment , which was 
equipped with the same superconducting magnet and 
calorimeter but with a different tracking system and with 
a gas RICH. For ascent and ffoat muons this systematic 
uncertainty decreased from ~ 9% at 0.3 GeV/c to about 
7% above 1 GeV/c. These uncertainties were not in- 
cluded in the data presented in the tables and the figures. 



IV. RESULTS 

We selected 37864 /z" and 47043 ^l+ between 0.2 and 
120 GeV/c at ground (1000 g/cm^); 5081 /i" between 



0.3 and 40 GeV/c and 2715 between 0.3 and 2 GeV/c 
during the ascent (7-850 g/cm^); 1601 /z~ between 0.18 
and 20 GeV/c and 2063 ii+ between 0.18 and 2 GeV/c 
at fioat altitude (3.3-4.6 g/cm^, mean atmospheric depth 
of 3.9 g/cm^). From these we obtained the muon fluxes 
(Jp- p+) according to: 



J„- 



1 



Tlive {X) X Gp- _ ^ 
xA^u-, u+(P,^), 



X e(P, X) X AP 



(1) 



where X is the atmospheric depth, Tuve is the live time, 



G,, 



are the geometrical factors for fjL and 



is the combined selection efficiency, AP is the width of 
the momentum bin corrected for ionization losses to the 
top of the payload, P the momentum and N^-^ ^+{P) 
is the selected number of /i^ and /i+ . The fractional 
live time decreased from 0.97240 ± 0.00001 at ground to 



0.2690 ±0.0006 at float altitude as indicated in Table 

Figure || and Table IV show the muon spectra at 
float, corresponding to 3.9 g/cm^ of residual atmosphere. 
These muon fluxes are interesting since these muons are 
the products of the first interaction between the primary 
cosmic rays and the air nuclei. Hence, along with the si- 
multaneous measurement of the primary spectra of pro- 
ton and helium nuclei ||^ these data provide a useful test- 
bench for studying the pion production in nucleon-air in- 
teraction used in the calculation of atmospheric showers. 
As discussed in section 



HIB2 



we limit the positive muon 
data to momenta below 2 GeV/c. The average muon 
charge ratio on this momentum interval is 1.59 ± 0.06. 

In Table IH we present the measured muon fiuxes at 
several atmospheric depths and momenta interval. The 
symbol FAD stands for Flux-weighted Average Depth 
p3[ obtained according to: 



FAD(P) 



jX{t)eu,e{t)J{P)dt 
Jeu,e{t)J{P)dt ' 



(2) 



where eu^e is the fractional live time. The depth and 
momentum intervals were chosen to match the pub- 
lished data of /i~ fiux growth curves by the MASS89 and 
MASS91 experiment [[o|j2|]. 

Figure |^ shows the flux growth curves for (a) negative 
and (b) positive muons for several momentum bins. For 
each momentum interval we fitted the data at large at- 
mospheric depths {X > 190 g/cm^) with an exponential 
function [ESl: 



J(P,X) = k{P)e 



-x/A(P) 



(3) 



where k and A are obtained from the fits. The resulting 
best fits are shown in Figure |^ as solid lines. As found 
in the MASS89 and MASS91 experiments a 
nearly linear relation exists between the attenuation 
length (A) and the mean momentum (P) in unit of GeV/c 
in the 190 to 1000 g/cm^ range. The relation resulting 
from CAPRICE94 measurements of /i^ is: 
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A[g/cm2] = (263 ± 14) + (150 ± 15) x P, 



(4) 



This relation holds also for the CAPRICE94 /i+ flux 
growth curves. It is worth pointing out that in deter- 
mining relation Q we used also the /i^ fluxes at ground. 
Equation |^ can be compared with the one determined 
from MASS89 fj,' data pi as: 



A[g/cm2] = (283 ± 24) + (93 ± 16) x P. 



(5) 



Both the above expressions reproduce the data within 
errors over the momentum range of these experiments, 
but will differ when extended to much larger momenta. 

Figure |l^ shows the relative difference between the 
/i~ fluxes obtained in this analysis and the MASS89 |2^] 
and MASS91 [0 experiments as a function of atmo- 
spheric depth. The comparison is done for muon mo- 
menta below 1 GeV/c (a) and between 1 and 2 GeV/c 
(b) . The dashed lines indicate the average difference be- 
tween this analysis and MASS89 and the solid lines the 
average difference with MASS91. Considering the errors 
in the data points, a good agreement is found in the 1 to 
2 GeV/c interval among different measurements. How- 
ever, below 1 GeV/c the CAPRICE94 results are signifl- 
cantly higher than the results from the MASS91 and, to a 
lesser extent, from the MASS89 experiments. These dif- 
ferences could be caused by solar activity or geomagnetic 
effects. In fact, the MASS89 experiment was launched 
from Prince Albert, Saschatcewan, Canada, during a 
period of maximum solar activity while MASS91 was 
launched from Fort Sumner, New Mexico, and flew at 
an average geomagnetic cutoff of about 4.5 GV/c. 

Geomagnetic effects also are observed in the muon 
charge ratio. Figure |l^ shows the /i"*" to fj,~ ratio as a 
function of atmospheric depth measured in the momen- 
tum intervals 0.3-1 GeV/c (a) and 1-2 GeV/c (b) by this 
experiment, by the recent CAPRICE98 experiment [ p5| , 
which flew from Fort Sumner on 28-29 May 1998, in the 
range 0.3-0.9 GeV/c by the HEAT95 experiment [|6), 
which flew from Lynn Lake on the 23rd of August 1995, 
and in the range 0.3-0.9 GeV/c (a) and 0.9-1.5 GeV/c 
(b) by the MASS91 experiment [Q. It can be seen that 
the CAPRICE94 low momenta charge ratios are higher 
than the MASS91 and CAPRICE98 ones. Moreover, the 
CAPRICE94 data show a dependence on the atmospheric 
depth, which is also visible in the HEAT data. Similar 
latitude effects also can be seen in Figure which shows 
(a) the CAPRICE94 data at float along with the charge 
ratios measured by the CAPRICE98 experiment and 
the MASS91 experiment |^8| and (b) the ground muon 
data reported here and from the CAPRICE97 experiment 
| p2| , which was carried out in Fort Sumner during Spring 
1997. 

Figure |l^ shows the measured spectra of negative 
muons for nine depth intervals. Above 1.5 GeV/c the 
spectra between 3.9 and 250 g/cm^(at larger atmo- 
spheric depths unacceptable power law flts were found 
for this momentum range) are power law in momentum 



with a fairly constant spectral index of —2.30 ± 0.04 
that can be compared with —2.5 ± 0.2 between 20 and 
400 g/cm2 above 2 GeV/c from MASS89 |2| and with 
-2.45 ±0.05 between 25 and 250 g/cm^ above 1.5 GeV/c 
from MASS91 ITotl. 



V. CONCLUSIONS 

In this paper we have presented new results on at- 
mospheric data measured with the CAPRICE94 exper- 
iment both for positive and negative muons. The data 
cover a large atmospheric depth range from close to the 
top of the atmosphere (3.9 g/cm^) down to ground level 
(1000 g/cm^). 

The data were compared with other experimental re- 
sults 1 10 2^ that were obtained using the same supercon- 
ducting magnet but with different identifying detectors. 
The muon spectra measured by the different experiments 
at high momenta (above 1 GeV/c) are in good agreement 
considering the overall uncertainty of the measurements 
(^ 10 — 15%). At lower energy, the comparison between 
the results of CAPRICE94 and those of MASS89/91 (the 
CAPRICE94 fi- fluxes are about 10-20% higher than the 
ones measured in the other two experiments) indicates 
solar modulation and geomagnetic effects. It is worth 
pointing out that the differences between the different 
measurements cannot account for the discrepancies found 
at low momenta, while comparing the experimental data 
with the theoretical calculation, which are in some cases 
as large as 70% (see pT[). 
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TABLE I. Selection used for identifying muons at different altitudes. See the text for a detailed discussion about each 
selection. 





Selection 


Ground 


Ascent 


Float 




{X = 1000 g/cm^) 


(3.9 <X < 1000 g/cm^) 


(X = 3.9 g/cm^) 


Tracking 


Used 


Used 


Used 


dE/dX 


Used 


Used 


Used 


/3tof 


Not used 


Not used 


Used 


Calorimeter 


Used for 0.3 < R < 3.2 GV/c 


Used for R > 0.3 GV/c 


Used for R > 0.3 GV/c 


RICH 


All criteria for 


All criteria for R < 0.53 GV/c, 


All criteria 




R < 0.53 GV/c 


and criterion 1 for 








R > 0.53 GY/c,X < 70 g/cm^ 






Selection 


Ground 


Ascent 


Float 




{X = 1000 g/cm^) 


(3.9 < X < 1000 g/cm^) 


(X = 3.9 g/cm^) 


Tracking 


Used 


Used 


Used 


dE/dX 


Used 


Used 


Used 


/3tof 


Not used 


Used 


Used 


Calorimeter 


Used for 0.3 < i? < 3.2 GV/c 


Used for R > 0.3 GV/c 


Used for R > 0.3 GV/c 


RICH 


All criteria 


All criteria 


All criteria 
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TABLE II. Sources for muon background 



Source 



Rejection criteria 



Residual contamination 



Albedo 
e~, e"*" 

Atmospheric mesons 

Locally produced 
mesons 



Protons^ 



Heavier than p nuclei^] 



time-of-fiight 
RICH and calorimeter 
RICH 

RICH, track 
reconstruction 



Aof and RICH; 
estimated and subtracted 
at X < 5 g/cm^ for 7? > 1.5 GV/c, 
at X = 1000 g/cm^ for > 3 GV/c 



dE/dX and RICH 



none 
negligible 
negligible for R < 0.4 GV/c, 
estimated < 2% for R > 0.4 GV/c 
at X < 5 g/cm^ 
< 20% for 0.5 < i? < 1 GV/c, 

< 10% for i? > 1 GV/c; 

at X > 7 g/cm^ 
negligible for R < 0.4 GV/c, 
< 1% 7£ > 0.4 GV/c 

at X < 5 g/cm^ 
negligible for R < 1.5 GV/c, 

< 1% for R > 1.5 GV/c; 

at 7 < X < 890 g/cm^ negligible; 
at X = 1000 g/cm^ 
negligible for R > 3 GV/c, 
< 1% for i? > 3 GV/c 
negligible 



^Only for p 
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TABLE III. Measured atmospheric growth fluxes (in units of particles/(GeV cm^ sr s)) for negative and positive muons 
in the 0.3 to 40 GeV/c momentum range. Results are given for the following momentum intervals: (I) 0.3-0.53 GeV/c, (II) 
0.53-0.75 GoV/c, (III) 0.75-0.97 GcV/c, (IV) 0.97-1.23 GoV/c, (V) 1.23-1.55 GcV/c, (VI) 1.55-2 GeV/c, (VII) 2-3.2 GeV/c, 
(VIII) 3.2-8 GeV/c, (IX) 8-40 GeV/c. The errors include both statistical and systematic errors. 



Depth interval 


A 






B 


C 




D 




Duration (s) 


85800 






610 


620 




560 






0.972 






0.949 


0.831 




0.646 




Initial depth (g/cm^) 


1000 






890 


580 




380 




Final depth (g/cm^) 


1000 






580 


380 




250 




I fi~ Flux 


(1.45 ± 0.04) X 


10" 


-y 


(3.49"^^,"^) X 10"^ 


(6.4 ± 1.3) X 10-^ 




(1.77 ±0.32) X 10" 


-'2 


FAD (g/cm^) 


1000.0 






709.0 


463.7 




308.0 




II Flux 


(1.36 ±0.03) X 


10" 


-3 


(2.55 ± 0.43) X 10"^ 


(6.67 ± 0.75) X 10" 


3 


(1.05 ±0.16) X 10" 


-2 


FAD (g/cm^) 


1000.0 






700.2 


466.4 




307.8 




III IJ,~ Flux 


(1.23 ± 0.02) X 


10" 


-3 


(1.94 ± 0.37) X 10"3 


(4.50 ± 0.60) X 10" 


3 


(7.26 ± 1.21) X 10" 


-3 


FAD (g/cm^) 


1000.0 






703.2 


467.0 




306.2 




IV n~ Flux 


(1.13 ±0.02) X 


10" 


-3 


(1.73 ±0.32) X 10"3 


(4.45 ± 0.55) X 10" 


3 


(5.35 ± 0.92) X 10" 


-3 


FAD (g/cm^) 


1000.0 






700.7 


469.0 




309.1 




Y fi- Flux 


(9.42 ± 0.17) X 


10" 


-4 


(1.34 ± 0.25) X 10"* 


(3.32 ±0.42) X 10" 


3 


(3.65 ±0.65) X 10" 


-3 


FAD (g/cm^) 


1000.0 






701.7 


470.1 




310.2 




VI /i" Flux 


(7.56 ±0.13) X 


10" 


-4 


(1.21 ±0.20) X 10"* 


(2.36 ± 0.30) X 10" 


3 


(2.59 ± 0.46) X 10" 


-3 


FAD (g/cm^) 


1000.0 






707.6 


471.8 




309.0 




VII fi- Flux 


(5.11 ±0.07) X 


10" 


-4 


(6.63 ± 0.92) X 10"" 


(1.27 ±0.14) X 10" 


3 


(1.89 ± 0.29) X 10" 


-3 


FAD (g/cm^) 


1000.0 






708.0 


469.2 




310.4 




VIII ^i" Flux 


(1.75 ± 0.02) X 


10" 


-4 


(2.32 ± 0.28) X 10-* 


(2.81 ± 0.32) X 10" 


4 


(4.34 ± 0.67) X 10" 


-4 


FAD (g/cm^) 


1000.0 






720.0 


472.3 




310.2 




IX ij,- Flux 


(1.33 ± 0.02) X 


10" 


-6 


(1.38 ±0.26) X 10"^ 


(1.97 ±0.33) X 10" 


5 


(1.37t°;**) X 10- 




FAD (g/cm^) 


1000.0 






715.6 


478.1 




309.1 




I /i+ Flux 


(1.76 ± 0.05) X 


10" 


-3 


(2.30+i: ?^,) X 10-* 


(8.5 ± 1.6) X 10-' 




(2.28 ± 0.40) X 10" 


-'2 


FAD (g/cm^) 


1000.0 






691.9 


459.2 




307.3 




II IJ.+ Flux 


(1.64 ±0.03) X 


10" 


-3 


(2.45l^:g) X 10-* 


(9.6 ± 1.4) X 10-* 




(1.70 ± 0.29) X 10" 


-2 


FAD (g/cm^) 


1000.0 






690.7 


462.6 




308.6 




III fi+ Flux 


(1.53 ±0.03) X 


10" 


-3 


(2.08li;:I*) X 10-* 


(4.85 ± 0.92) X 10" 


3 


(1.29 ± 0.23) X 10" 


-2 


FAD (g/cm^) 


1000.0 






703.2 


462.4 




307.3 




IV /i+ Flux 


(1.31 ±0.02) X 


10" 


-3 


{1.83t°o ll) X 10-* 


(5.49 ± 0.89) X 10" 


3 


(1.05 ±0.19) X 10" 


-2 


FAD (g/cm^) 


1000.0 






696.9 


463.7 




309.2 




V Flux 


(1.12 ±0.02) X 


10" 


-3 


(9.7tt.?,) X 10-* 


(3.49 ±0.63) X 10" 


3 


(5.59 ±1.13) X 10" 


-3 


FAD (g/cm^) 


1000.0 






692.6 


464.1 




309.3 




VI /i+ Flux 


(9.36 ± 0.16) X 


10" 


-4 


(1.23l°S) X 10-* 


(2.00 ± 0.40) X 10" 


3 


(3.84 ± 0.78) X 10" 


-3 


FAD (g/cm^) 


1000.0 






712.4 


468.1 




310.7 





Depth interval 


E 




F 


G 


H 


Duration 


740 




640 


660 


690 




0.516 




0.449 


0.396 


0.358 


Initial depth (g/cm^) 


250 




190 


150 


120 


Final depth (g/cm^) 


190 




150 


120 


90 


I IJ,~ Flux 


(1.44 ±0.19) X 10" 


-2 


(2.06 ± 0.27) X 10-^ 


(1.70 ± 0.25) X lO-'' 


(1.69 ± 0.25) X lO-'' 


FAD (g/cm^) 


218.8 




173.8 


135.3 


104.1 


II ^- Flux 


(1.26 ±0.11) X 10" 


-2 


(1.26 ±0.13) X 10-^ 


(1.21 ±0.13) X 10-2 


(1.21 ±0.13) X 10-2 


FAD (g/cm^) 


219.1 




174.1 


135.1 


104.3 


III fi- Flux 


(1.05 ±0.10) X 10" 


-2 


(8.6 ± 1.0) X 10-* 


(9.1 ± 1.1) X 10-* 


(8.7 ± 1.1) X 10-* 


FAD (g/cm^) 


219.3 




174.3 


135.0 


104.1 


IV Flux 


(6.51 ±0.73) X 10" 


-3 


(6.35 ±0.83) X 10"* 


(6.31 ±0.86) X 10"* 


(6.05 ±0.86) X 10-* 


FAD (g/cm^) 


219.1 




174.1 


135.1 


104.2 


V fi- Flux 


(4.03 ± 0.52) X 10" 


-3 


(6.03 ± 0.73) X 10-* 


(5.15 ±0.71) X 10-* 


(4.55 ± 0.67) X 10-* 


FAD (g/cm^) 


218.3 




173.7 


135.5 


104.4 


VI fi- Flux 


(3.36 ± 0.40) X 10" 


-3 


(3.42 ± 0.46) X 10-* 


(3.70 ± 0.51) X 10-* 


(2.92 ± 0.46) X 10-* 


FAD (g/cm^) 


219.0 




173.9 


135.3 


104.5 


VII fi- Flux 


(1.65 ±0.17) X 10" 


-3 


(1.49 ± 0.19) X 10-* 


(1.85 ± 0.22) X 10-* 


(1.19 ±0.18) X 10-* 


FAD (g/cm^) 


219.8 




173.8 


135.4 


104.6 
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VIII Flux 


(3.78 ± 0.42) X 10" 


-4 


(3.66 ± 0.47) X 10 " 


(3.51 ± 0.49) X 10 * 


(3.63 ± 0.50) X 10 


FAD (g/cm^) 


219.6 




174.2 


135.1 


104.4 


IX Flux 


(2.43 ± 0.41) X 10" 


_5 


(i-38i;;:^?) X 10-' 


(1.42l»:«) X 10 5 


{1.89toil) X 10 5 


FAD (g/cm^) 


219.8 




174.9 


134.6 


104.6 


I /i+ Flux 


(2.07 ±0.25) X 10" 





(2.57 ±0.32) X lO""" 


(2.09 ± 0.30) X lO"'' 


(1.85 ± 0.28) X 10" 


FAD (g/cm^) 


219.0 




174.0 


135.5 


104.0 


II n+ Flux 


(1.68 ±0.18) X 10" 




(1.96 ±0.23) X 10"^ 


(1.52 ±0.21) X 10"^ 


(2.06 ±0.25) X 10" 


FAD (g/cm^) 


219.0 




174.2 


135.0 


104.0 


III n+ Flux 


(1.17 ±0.15) X 10" 




(1.39 ±0.19) X 10"^ 


(1.71 ± 0.22) X 10"^ 


(1.24 ± 0.19) X 10" 


FAD (g/cm^) 


219.1 




173.4 


135.2 


104.6 


IV Flux 


(9.1 ± 1.2) X 10"' 




(9.3 ± 1.4) X 10"^ 


(1.05 ±0.15) X 10"^ 


(9.0 ± 1.5) X 10"^ 


FAD (g/cm^) 


219.5 




173.8 


135.1 


104.4 


V Flux 


(5.44 ± 0.83) X 10" 


-3 


(7.4 ± 1.1) X 10"^ 


(8.0 ± 1.2) X 10"^ 


(5.3 ± 1.0) X 10"^ 


FAD (g/cm^) 


218.7 




173.4 


135.5 


104.3 


VI /i+ Flux 


(2.68 ± 0.49) X 10" 


-3 


(3.56 ±0.65) X 10"^ 


(3.67 ±0.69) X 10"'^ 


(4.05 ± 0.74) X 10" 


FAD (g/cm^) 


219.2 




173.5 


134.8 


104.2 



Depth interval 


I 




J 




K 




L 




Duration 


720 




1300 




1660 




60520 




eiive 


0.327 




0.303 




0.289 




0.269 




Initial depth (g/cm^) 


90 




65 




35 




3.3 




Final depth (g/cm^) 


65 




35 




15 




4.6 




I n' Flux 


(1.66 ±0.26) X 10" 


•2 


(1.01 ±0.14) X 10" 


-2 


(4.89 ±0.86) X 10" 


-:-! 


(1.40 ±0.11) X 10" 


-3 


FAD (g/cm^) 


77.3 




50.7 




25.7 




3.9 




II /i" Flux 


(1.01 ±0.12) X 10" 


2 


(8.58 ±0.89) X 10" 


-3 


(5.56 ± 0.65) X 10" 


-3 


(1.04 ± 0.06) X 10" 


-3 


FAD (g/cm^) 


77.2 




49.8 




25.1 




3.9 




III Flux 


(8.9 ± 1.2) X 10"^ 




(5.81 ± 0.73) X 10" 


-3 


(3.69 ± 0.52) X 10- 


-3 


(6.76 ± 0.51) X 10" 


-4 


FAD (g/cm^) 


77.2 




50.3 




25.1 




3.9 




IV /i" Flux 


(5.59 ± 0.84) X 10" 


3 


(3.92 ± 0.55) X 10" 


-3 


(1.73 ± 0.33) X 10- 


-3 


(4.68 ± 0.38) X 10" 


-4 


FAD (g/cm^) 


77.2 




50.6 




25.9 




3.9 




V /i" Flux 


(3.65 ± 0.61) X 10" 


3 


(2.64 ± 0.41) X 10" 


-3 


(1.61 ±0.29) X 10" 


-3 


(2.76 ± 0.26) X 10" 


-4 


FAD (g/cm^) 


77.4 




50.1 




25.2 




3.9 




VI M" Flux 


(2.41 ± 0.42) X 10" 


3 


(2.08 ± 0.31) X 10" 


-3 


(7.6 ± 1.7) X 10" 


i 


(1.93 ±0.19) X 10" 


-4 


FAD (g/cm^) 


77.2 




50.4 




26.2 




3.9 




VII ^i" Flux 


(1.17 ±0.18) X 10" 


3 


(7.4 ± 1.1) X 10"' 


1 


(4.98 ± 0.84) X 10- 


-4 


(9.19 ± 0.79) X 10" 


-5 


FAD (g/cm^) 


77.2 




50.3 




25.0 




3.9 




VIII Flux 


(2.61 ± 0.43) X 10" 


4 


(1.58 ± 0.26) X 10" 


-4 


(1.16 ±0.20) X 10" 


-4 


(1.76 ± 0.17) X 10" 


-5 


FAD (g/cm^) 


77.6 




50.3 




24.9 




3.9 




IX /X" Flux 


(7.7 ±4.2) X 10"® 


(imt'ofe) X 10- 


5 


(4.4t?:^) X 10"8 




(1.70 ± 0.20) X 10" 


-6 


FAD (g/cm^) 


77.7 




49.5 




26.0 




3.9 




I Flux 


(2.34 ±0.33) X 10" 


■2 


(1.60 ±0.19) X 10" 


-2 


(6.4 ± 1.0) X 10" 




(2.32 ±0.18) X 10" 


-3 


FAD (g/cm^) 


76.9 




50.7 




26.1 




3.9 




II /i+ Flux 


(1.53 ±0.22) X 10" 


2 


(1.00 ±0.12) X 10" 


-2 


(6.63 ± 0.91) X 10- 


-3 


(1.83 ±0.11) X 10" 


-3 


FAD (g/cm^) 


77.5 




50.2 




25.1 




3.9 




III Flux 


(1.07 ±0.18) X 10" 


2 


(5.08 ± 0.86) X 10" 


-3 


(4.05 ± 0.69) X 10- 


-3 


(1.08 ± 0.08) X 10" 


-3 


FAD (g/cm^) 


77.6 




50.6 




24.7 




3.9 




IV At+ Flux 


(7.8 ± 1.4) X 10"^ 




(5.10 ±0.79) X 10" 


-3 


(3.80 ± 0.62) X 10" 


-3 


(5.96 ± 0.49) X 10" 


-4 


FAD (g/cm^) 


77.4 




50.1 




24.8 




3.9 




V /x"^ Flux 


(6.3 ± 1.1) X 10"^ 




(2.50 ± 0.49) X 10" 


-3 


(2.29 ± 0.43) X 10" 


-3 


(4.49 ± 0.38) X 10" 


-4 


FAD (g/cm^) 


77.3 




50.8 




24.4 




3.9 




VI Ai+ Flux 


(3.15 ± 0.66) X 10" 


3 


(1.99 ±0.37) X 10" 


-3 


(1.07l«;^i) X 10" 


3 


(2.91 ±0.28) X 10" 


-4 


FAD (g/cm^) 


77.5 




50.5 




25.5 




3.9 
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TABLE IV. The /i^ and fx^ fluxes at float altitude (3.9 g/cm^ of mean residual atmosphere). Columns 1 and 2 are, 
respectively, the rigidity bins in the spectrometer and the mean momenta at the top of the payload. Columns 3 and 4 give the 
resulting fj,~ and /i"*" fluxes, column 5 the /i"*" to fj,~ ratio. The errors include both statistical and systematic errors. 



Rigidity 
bin 
GV/c 


Moan 
momentum 
GeV/c 


Flux 
(GeV cm^ sr 




M"*" / M 


0.15 - 0.20 


0.19 


(1.66+^8) xlO-" 


(2.88 ± 0.80) xlO"^ 


i-74l"o;?^ 


0.20 - 0.30 


0.27 


(1.54 ± 0.20) xlO"^ 


(3.05 ± 0.35) xlO"* 


1.97 ± 0.35 


0.30 - 0.40 


0.37 


(1.56 ± 0.18) xlO"^ 


(2.50 ± 0.23) xlO"^ 


1.60 ± 0.24 


0.40 - 0.60 


0.51 


(1.35 ± 0.10) xlO~^ 


(2.33 ± 0.12) xlO"^ 


1.73 ± 0.16 


0.60 - 0.75 


0.69 


(8.97 ± 0.76) xlO""* 


(1.65 ± 0.10) xlO"^ 


1.84 ± 0.19 


0.75 - 0.90 


0.84 


(7.64 ± 0.70) xlO"* 


(1.11 ± 0.08) xlO"^ 


1.45 ± 0.17 


0.90 - 1.10 


1.01 


(5.23 ± 0.49) xlO"* 


(7.18 ± 0.52) xlO"* 


1.37 ± 0.16 


1.10 - 1.30 


1.21 


(3.94 ± 0.42) xlO~* 


(5.78 ± 0.47) xlO"'' 


1.47 ± 0.20 


1.30 - 1.60 


1.46 


(2.43 ± 0.27) xlO"* 


(3.83 ± 0.31) xlO"'' 


1.58 ± 0.21 


1.60 - 2.0 


1.80 


(1.73 ± 0.19) xlO"^ 


(2.75 ± 0.25) xlO"* 


1.60 ± 0.23 


2.0 - 3.0 


2.46 


(9.87 ± 0.96) xlO"^ 






3.0 - 5.0 


3.87 


(3.37 ± 0.39) xlO"^ 






5.0 - 10.0 


6.97 


(8.60 ± 1.22) xlO"'' 






10.0 - 20.0 


13.90 


(2.07 ± 0.42) xlO"*^ 
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FIG. 1. Schematic view of the CAPRICE94 apparatus. 
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FIG. 2. Cosine zenith angle distribution normalized to the total number of events for rimons of both signs selected between 
0.15 and 2 GV/c at ground level (dotted histogram), during the ascent (dashed histogram) and at float (solid histogram). 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Rigidity (GV/c) 

FIG. 3. The distribution of 13 from the timc-of-fiight information as a function of rigidity for float data. The solid line is 
the lower limit of the muon selection. The distributions corresponding to the various particles are labelled according to the 
particle's species. The figure comprises about 379000 events of which about 347000 are protons, 9000 deuterons and 6000 
helium nuclei. 
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FIG. 4. Total detected energy loss in the calorimeter divided by the rigidity of all particles versus rigidity for float data. A 
negative sign is assigned to the rigidity of negative particles. The two dense band are due to non-interacting particles. The 
solid line at 60 indicates the chosen upper limit used to select muons. The figure comprises about 615000 events of which about 
16000 are negative particles. 
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FIG. 5. Measured RICH Cherenkov angle as a function of rigidity for float data. A negative sign is assigned to the rigidity 
of negative paxticles. The figure comprises about 291000 events of which about 8000 are negative particles. The solid lines 
include the events accepted as muons in the RICH selection. 
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FIG. 6. Geometrical factor as a function of deflection for ground and flight analysis. The lines indicate one standard deviation 
confidence interval. 
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FIG. 7. Muon selection efficiencies for float data as a function of rigidity. The hatched area indicates one standard deviation 
confidence interval of the combined efficiency. 
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FIG. 8. The (full circles) and /i (open circles) fluxes at float (3.9 g/cm of mean residual atmosphere) as a function of 
momentum. 
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FIG. 9. Atmospheric growth curves for a) fi in the left panel and b) /i^ in the right panel. From top to bottom are the 
momentum ranges in GeV/c: 0.3-0.53 (scaled by lO''), 0.53-0.75 (10*), 0.75-0.97 (10^), 0.97-1.23 (10^), 1.23-1.55 (10), 1.55-2 
(1), 2-3.2 (1), 3.2-8 (1) and 8-40 (1). The fi'^ results are shown up to 2 GeV/c. The solid lines are exponential fits according 
to equation H. 
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FIG. 10. Relative difference between the ii~ fluxes obtained in this analysis and the MASS89 |23| and MASS91 experi- 
ments as a function of atmospheric depth. The comparison is done for muon momenta below 1 GeV/c in the upper panel (a) 
and between 1 and 2 GeV/c in the lower panel (b). The average difference between CAPRICE94 and MASS89 is shown as a 
dashed line and between CAPRICE94 and MASS91 as a solid hue. 
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FIG. 11. Muon charge ratio as a function of atmospheric depths in two momentum intervals: (a) 0.3-1 GeV/c in the upper 
panel and 1-2 GeV/c in the lower panel (b), obtained in this analysis and from the CAPRICE98 [|||, HEAT95 ||§ and MASS91 
experiments |l^; the latter data refer to the momentum bins 0.3-0.9 GeV/c and 0.9-1.5 GeV/c. 
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FIG. 12. Muon charge ratios as a function of momentum. In the upper panel (a) are shown the ratios obtained in this 
analysis at 3.9 g/cm'^and those measured by the CAPRICE98 experiment p7| at 5.5 g/cm'^and with the MASS91 experiment 
[ psi at 5.7 g/cm^. In the lower panel (b) are shown the ratios measured at ground by this experiment (at 1000 g/cm^) and by 
the CAPRICE97 experiment (at 886 g/cm^) 



27 



10' 



CO : 

(0 1 ^ 



> 

O 10 



CO . 

o 



O 10 

X _. 



xlO' • 

x10' 



o 10^^ x10' 



^ 10 - 



10 - 



-4 

10 - 



o 



x10' □ 

X10' A 
X10 A 

x1 ▼ 



o 



o 



□ 



A 



T 





• 


• 




o 


o 


• 








■ 


■ 


■ 



□ 



A 



T 





□ 



□ 



A ▲ 



A 



* * * ^ 




-5 " 

10 - 



10 



Momentum (GeV/c) 



FIG. 13. Negative muon spectra for several atmospheric depth intervals. From top to bottom are the atmospheric depth 
ranges in g/cm^: 3.3-4.6 (scaled by 10'^), 7-25 (scaled by 10^), 25-70 (10^), 70-115 (10^), 115-165 (10^), 165-250 (10), 250-350 
(1), 350-850 (1) and 1000 (1). The solid lines are power law fits. 
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